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Abstract. Collective effects, such as waves and instabilities, are integral to our 
understanding of most plasma phenomena. We have been able to study these in 
ultracold neutral plasmas by shaping the initial density distribution through spatial 
modulation of the ionizing laser intensity. We describe a relay imaging system for 
the photoionization beam that allows us to create higher resolution features and its 
application to extend the observation of ion acoustic waves to shorter wavelengths. 
We also describe the formation of sculpted density profiles to create fast expansion of 
plasma into vacuum and streaming plasmas. 
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1. Introduction 

Collective modes are central to our understanding of plasma dynamics and transport 
properties [TJ [2]. They often lead to instabilities [3] that can dominate plasma 
dynamics and create significant complications for practical applications such as magnetic 
confinement fusion. In principle, and for most regimes of plasma density and 
temperature, these phenomena are well described within classic theoretical frameworks 
such as continuum fluid or kinetic descriptions, but the possible complexity in plasmas 
is continually leading to the observation of interesting new phenomena. In addition, 
at relatively high densities and/or low temperatures, when plasmas become strongly 
coupled, or non-ideal [I], plasma dynamics becomes hard to describe theoretically or 
simulate numerically, so new experiments probing this regime are of great fundamental 
interest and important for understanding the behavior of matter under extreme 
conditions of energy and density [5]. 

Recent advances in plasma creation have allowed for the excitation of ion collective 
modes in ultracold neutral plasmas (UCNPs) [6]. UCNPs [TJ [8] are formed by 
photoionizing laser cooled atoms near the ionization threshold, and they access an exotic 
regime of plasma physics in which electron and ion temperatures are orders of magnitude 
colder than in traditional neutral plasmas. These plasmas have unique properties arising 
from their ultracold temperature, small size, and expansion dynamics. They also display 
correlated particle dynamics [9j [101 HU [121 1131 E] that reflect strong coupling [2] of the 
constituents, which allows them to probe a regime of nature in which correlated, many- 
body interactions are dominant. 

In strongly coupled systems, the Coulomb interaction energy exceeds the thermal 
energy, as characterized by a Coulomb coupling parameter, T > 1, where 



Here T is the temperature and a is the Wigner-Seitz radius, a = [3/ (47m)] 1 / 3 for density 
n. Because of the low temperature in UCNPs, ions equilibrate with 2 < Tj < 5 [TU| H~2]. 
It is possible to set initial conditions for large electron T e , but heating processes involving 
many different collisional effects will lower the coupling to T e < 0.2 [TB"| ITo^ [T7J [TS1 US] 
within a few microseconds. For T > 1, particle correlation effects become important. 

A wide range of phenomena has been studied in UCNPs, such as expansion of the 
plasma into surrounding vacuum with [20] and without [2TJ [HJ [22] an applied magnetic 
field, adiabatic cooling of electrons during expansion [20l [T9] , three-body recombination 
at ultracold temperatures [231 H21 [21] , spontaneous evolution of a frozen Rydberg gas 
into an ultracold plasma [22], equilibration of ions after plasma creation [TT| 126] . and 
disorder induced heating of electrons [HI [T7] and ions [TTJ [2S1 El EH] and kinetic energy 
oscillations [HI [26J [2Z] due to effects of strong coupling. Langmuir oscillations have 
been resonantly excited and used as a probe of plasma density [28, 29J, and an E x B 
electron drift instability was also observed [30] . 

Recently, Castro et al. \6\ used a new technique involving spatially masking the 
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photoionization-laser intensity to shape the initial density distribution of UCNPs and 
excite long- wavelength ion acoustic waves (IAWs). The excitation wavelength was 
as short as A ~ 500 with the limit set by diffraction effects. Shorter wavelength 
excitations would access more interesting physics, such as the transition from a sound 
wave to a dispersionless ion plasma oscillation when the wavelength becomes on the 
order of the Debye screening length, X D = (eo/ceTe/nee 2 ) 1 / 2 , where T e and n e are 
the electron temperature and density respectively. Furthermore, corrections to the 
IAW dispersion relation due to strong coupling of ions are predicted for A < a 
[3T| |32| |33| |34"| |35| |36| [37] and sharp features in the ion density profile should also 
lead to large ion accelerations that may excite shock waves [H2 [7J and solitons [38J. 
In this paper we present improvements to the ionization system that decrease the 
observable excitation wavelength to ~100 //m, approaching the Debye length under 
accessible UCNP conditions and only 10 times the Wigner-Seitz radius. 

2. Experimental Setup 

Strontium UNCPs are created in a magneto-optical trap (MOT) containing ~ 300M 
laser cooled 88 Sr atoms by photoionizing approximately half of the atoms near the 
ionization threshold with a pulsed narrowband laser [7J. The density distribution of the 
atoms is a spherical gaussian, oc exp(— r 2 /2er 2 ), with characteristic radius a ~ 1 — 2 mm, 
and if the ionization beam intensity profile is uniform, the plasma has a similar form. 
The electrons have a tunable temperature, (^E e /ks ~ 1 — 100 K), determined by the 
energy difference (E e ) between the ionizing photon and the ionization potential. The 
ions, with their relatively large mass, inherit the temperature of the laser-cooled atoms 
(~ lOmK). 

Two-photon ionization is realized with one photon from a pulse amplified beam 
of light from the cooling laser, driving the Sr 1 So - x Pi transition at 461 nm, and a 
second photon from a pulsed dye laser, at 412 nm, which promotes the electron to 
the continuum. We use a mask on the second ionizing beam in order to spatially 
modulate the laser intensity. This simple technique allows one to spatially manipulate 
the ionization fraction, in turn controlling the initial density distribution of the UCNP. 
For instance masking half the beam with a straight edge produces a single Gaussian 
hemisphere of plasma. Further flexibility is provided by the option of reflecting the 
ionizing beam through the chamber a second time to increase the ionization fraction, 
and by masking only the second pass of the 412 nm ionizing beam we can create smaller 
amplitude perturbations. For example, a 100% contrast square wave mask was used 
in this fashion in previous studies to excite low amplitude IAWs with user-selected 
wavelength [6]. Although highly versatile, this straightforward technique is limited in 
terms of achievable feature size due to laser diffraction over the distance from the mask 
to the plasma. To overcome this limitation we have implemented a relay lens system to 
image the mask onto the plasma with much higher resolution. 

Our design consists of a matched, positive focal length lens pair (Thorlabs AC508- 
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Figure 1. Experimental setup showing the optical arrangement for ionizing and 
imaging an UCNP. The ionizing relay lenses are set at a distance for optimum imaging 
of the transmission mask onto the plasma, and the transmission masks are mounted 
in rotation and translation stages. To increase the ionization fraction the ionizing 
beam can be rerouted for a second pass through the chamber (not shown). Placing the 
mask in the beam path just before the second pass creates smaller amplitude density 
perturbations. 



200-A / = 200 mm achromatic doublets) in a 1:1 conjugate ratio. To achieve maximum 
performance, many factors, such as lens aberrations, alignment, and aperture must be 
balanced against experimental design limitations, but this configuration was predicted 
with Lambda Research's OSLO optical design software to produce a contrast transfer 
function (CTF) of 50% for a 50cyc/mm pattern. The lenses are held in place with a 
custom lens tube adapter that can be securely attached to the vacuum chamber. This 
direct mounting method ensures correct and reproducible mask-lens-plasma alignment 
(coaxial with ionizing beam and distance / from the plasma). Transmission masks are 
mounted a distance / back in a 6 Z rotation mount atop an xz translation stage (laser 
propagation k/\k\ = z). Rotational alignment ensures that ionized edges are aligned 
vertically with the plasma imaging system's line of sight, x alignment allows centering 
of a pattern transversely on the plasma, while micrometer z adjustment allows fine 
positioning of the mask image so that the atoms are centered within the depth of field. 

3. Improved Resolution 

We measured the contrast transfer function of the ionization relay system by imaging 
a 1951 USAF resolution target onto a CCD camera. The CTF is calculated from the 
normalized signal (S) as CTF(v) = (S max (u) - S min (u)) / (S max (u) + S min (u)), where v 
is the spatial frequency of the square wave pattern. The system exceeded design goals, 
with a CTF of ~ 100% at 20cyc/mm over the typical depth of field of the plasma 
±2 mm and maximum focal plane resolution of CTF ~ 100% at >80cyc/mm. Note, 
this later measurement was limited by CCD pixel size. 

In situ diagnostics of the ionizing relay system was accomplished with fluorescence 
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imaging of the Sr + ions, on the principal transition, 2 S;l/ 2 — 2 Pi/2 at 422 nm [39J (Fig. 
[1]). A sheet of light, with a beam waist of 625 /im, propagating perpendicular to the 
imaging system line of sight and 106° away from the 412 nm ionizing beam, is used 
to excite fluorescence. Relative plasma density in the illuminated region is measured 
directly by scanning the imaging laser over the doppler broadened resonance of the 
entire plasma cloud and summing the captured fluorescence in post analysis. Absolute 
calibration of density is accomplished with complimentary absorption imaging [TO] . 
From the 2D images we create ID slices by cropping a narrow strip perpendicular to 
the density features and integrating along the narrow dimension (Figs. 2(a) and 2(b)). 
For these measurements, a single pass of the ionizing beam through the chamber is 
used, which should lead to 100% density contrast for well-resolved patterns. Using these 
analysis tools and square wave masks of various wavelengths we measured contrast of the 
imaged plasma density, which represents a CTF of the entire ionizing/imaging system. 
Figure 2(b), shows one such measurement where we achieve our highest resolution of 
CTF=13.5% at lOcyc/mm. 

This resolution limit is significantly worse than what was measured for the 
photoionization relay during bench testing. Effects of plasma dynamics and 
misalignment of the planes of constant density with the line of sight of the imaging 
system cannot explain the discrepancy. We suspect that the imaging system is 
currently limiting our resolution and the plasma density features have higher contrast. 
Improvements to the imaging system are now underway. 



4. Experiments with Sculpted UCNPs 

The observed resolution represents a five-fold improvement over the results 
presented in [6], and we can now study plasmas with density modulations with spatial 
frequency as high as lOcyc/mm. To isolate the density modulations, density profiles 



such as the one shown in Fig. 2(b) are fit to a Gaussian profile and then divided by 



the fitted Gaussian. With time resolved imaging, we observed the evolution of the ion 



density modulation (Fig. 3(a) ), which can be described in terms of IAWs similar to [6]. 

IAWs are longitudinal density modulations with slow enough ion dynamics for electrons 

to remain isothermal. In an infinite homogeneous medium, IAWs have the following 

dispersion relation for frequency u and wave vector k = 2n/\, 

'lu\ 2 k B T e /M (2) 



kj l + ^Af,' 

where M is the ion mass. 

For the short- wavelength IAWs that are now accessible, the higher frequency allows 
a simpler analysis than used previously [6]. The oscillation period is small compared 
to the plasma expansion timescale, therefore the electron temperature and IAW 
wavelength, which affect the oscillation frequency, do not change significantly during 
the wave period. This permits a simple description of the density perturbation as a 
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Figure 2. Fluorescence images and density profiles, (left) On resonance 
fluorescence image of UCNP with lOcyc/mm modulation, 100 ns after plasma creation. 
Limitations to optical access of the ionizing beam into the vacuum chamber produced 
the resulting oblique orientation of the plasma, (right) Density profile of plasma (solid) 
with Guassian fit (dotted). The lack of contrast on the left side results from reflections 
of the ionizing beam inside the chamber. X 1 is the distance along a perpendicular to 
the density features. 



constant wavelength, constant frequency, standing IAW, i.e 5n = Ave' 1 * cos(kx)cos(ut) , 
where Aq is initial amplitude, and 7 is an exponential decay rate. A fit of this model to 
data shown in Fig. 3(a) yields a radial oscillation frequency of 00 — 3.6(1) x 10 6 s _1 and 
a decay rate 7 = 9(1) x 10 5 s -1 . 

Figure 3(b) shows the comparison of measurements and theory for the dispersion 
relation. At higher resolution, the wavevector is in the transition regime from linear, 
acoustic behavior to the limiting case of a dispersionless ion plasma oscillation with 



Upi = yne 2 /Me® for large k. With higher electron temperature, which is easily 
accessible, we can push further into this regime. The measurement, however, does 
not show the expected rollover. This could be due to the initial density modulation 
being much larger than in previous studies [6] because equation |2] is only valid when the 
excitation is a small perturbation to the density distribution. 

The observation of fast damping of the oscillation is also an interesting observation. 
The damping is much stronger than observed at longer wavelengths [6], which may 
indicate the influence of Landau damping since this effect becomes large as the 
wavelength decreases [40] . Variation of frequency across the sample due to variation 
of density and A# may also contribute to the decreasing amplitude. With our increased 
resolution, we can now study damping in detail. 



Beyond IAW dispersion and damping, we can also study a variety of situations in 
which plasmas are strongly accelerated into the vacuum or into another plasma. Large 
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Figure 3. IAW Evolution for lOcyc/mm perturbation and T e = 25 K. (left) 
Evolution of the normalized perturbation. Note the initial periodicity closely matches 
that of the 10 cyc/mm mask, with dashed lines indicating positions of amplitude peaks. 
Between 400 and 500 ns the nodes and anti- nodes invert, which marks the first quarter 
cycle for the standing wave. A more rigorous fit, as described in the text, yields a 
quarter period for the IAW of T/4 = 440 ± 40 ns. (right) Extended IAW dispersion 
plot. The solid line corresponds to the theory curve, with linear extension of large 
wavelength regime shown as a dotted line to emphasize the screening correction. Circles 
indicate previously measured data [6], while the asterisk shows the new measurement. 



plasma density gradients, accessible with high resolution photoionization, produce large 
plasma accelerations and velocities, enabling the study of nonlinear effects such as shock 
formation and solitons [31 [15], [3, EE]. The importance of the density gradient can be 
seen from a hydrodynamic description using the coupled ion-electron momentum balance 
equations at early times when ion velocities are small [28]. Acceleration of the ions (ao) 
is provided by the electron pressure (P e ) according to Mndo m — VP e ~ —V [nfcsT e (0)]. 
T e (0) is expected to be close to uniform [15], leading to acceleration that varies as Vn. 
For an unperturbed UCNP, the density profile is spherical gaussian with | — Vn| « n/a Q , 
giving the characteristic radial acceleration on a timescale of about 10 fis that has been 
studied extensively in the context of self-similar plasma expansion [7J. However, when 
sharp features with length scale / cr are imprinted on the plasma density, the local 
gradient and resulting plasma acceleration and velocity can be much larger. 

A simple configuration is the total attenuation of half the ionizing beam. This 
creates a single, hard-edged, hemisphere with rapid expansion perpendicular to the 
edge. Figure H] shows this configuration and a clear picture of the maximum density 
gradient currently observable. Further studies will allow us to characterize this effect 
and possibly observe strongly supersonic velocities that may produce shock waves. 

Another interesting configuration utilizes a single opaque wire for the transmission 
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Figure 4. Single hemisphere (left) 2D image of a plasma hemisphere. Half of the 
ionizing beam is blocked to produce one hemisphere of an otherwise spherical Gaussian 
plasma, (right) ID density profile of a plasma hemisphere. Preliminary studies show 
how the strong gradient produced here results in ions that are accelerated strongly 
into the vacuum. 




Figure 5. Streaming plasma (left) 2D image of two plasma hemispheres. An opaque 
wire is used as a transmission mask, producing a gap in the middle of the plasmas, 
(right) ID density profile of two plasma hemisphere. Ions are accelerated strongly into 
the void were they can interact and penetrate into the opposite hemisphere. 



mask, thus creating two plasma hemispheres with an adjustable gap in between (See 
Fig. [5]). In this geometry, ions accelerate into the central void, where they can interact. 
Preliminary studies have shown a strong dependance on experimental parameters 
(density, gradient and electron temperature) of the ability for streaming plasmas to 
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penetrate each other. Additional studies will probe subsonic and supersonic plasma 
velocities and look for effects of strong coupling. 

5. Conclusion 

A new relay lens system on the photoionization laser beam has resulted in at least a five- 
fold improvement in ionizing resolution compared to previous studies of sculpted UCNPs 
[6] , allowing access to many interesting plasma experiments. We have measured shorter 
wavelength IAWs with very fast damping and observed ions that are accelerated strongly 
by a sharp plasma gradient. Resolution of the imaging system may currently limit 
the observable resolution, suggesting the underlying density features are significantly 
sharper. 
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